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We study the electromagnetic force exerted on a pair of parallel slab waveguides by the light 
propagating through them. We have calculated the dependence of the force on the slab separation 
£ — ■ by means of the Maxwell-Stress tensor formalism and we have discussed its main features for the 

different propagation modes: spatially symmetric (antisymmetric) modes give rise to an attractive 
(repulsive) interaction. We have derived the asymptotic behaviors of the force at small and large 
separation and we have quantitatively estimated the mechanical deflection induced on a realistic 
( — 1 ' air-bridge structure. 
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I. INTRODUCTION 

CZ2 ' 

q ■ When two objects in close proximity are illuminated by a light source an optical force is exerted on each of them, 
whose sign can be either attractive or repulsive depending on the geometry of the objects and of the optical mode in 
Oh which light propagates. Physically, this force originates from the interaction of the induced dipoles in the dielectric 
media by the electromagnetic field of the light wave. Its magnitude is proportional to the light intensity and depends 
r/3 1 on the actual profile of the electromagnetic field and on its polarization. 

The recent advances of nanotechnologies have led to the realisation of solid-state samples whose sizes and separations 
, are so small that the optical forces can have a significant impact on the shape and the position of the object. In 
particular, one may expect in the next future to take advantage of these optical forces to control the growth and the 
self-assembling of artificial materials such as photonic crystals or random sphere assembly. Moreover optical forces 
may be also used to engineer the quantum state of the mechanical motion of nano-objects [l[ . 

In the last years, a significant amount of experimental and theoretical studies have concerned coupled resonant 
systems, such as spherical and disk-shaped whispering gallery resonators Q as well as double-layer photonic-crystal- 
slab cavities Q. Thanks to the extremely high quality factor of these systems, the light intensity in the resonator is 
in fact strongly enhanced as compared to the input power, which leads to a corresponding increase in the magnitude 
I of the force. 

, At the same time the case of two coupled air bridge silicon waveguides with square cross section has been investigated 
in from the reported numerical calculations, it turns out that the displacement of the silicon wire-bridges under the 
action of the induced force can reach values measurable with standard Atomic Force Microscope techniques already 
■ at reasonable input laser power. 

Motivated by this intense research and by the fast advances in the nanotechnological expertise in manipulating these 
systems, we here present a systematic characterization of the optical force between two parallel planar waveguides. 
For such a simple geometry most of the results can be obtained by analytical means, which provides useful insight 
k^~>! into the basic physics of the optical forces. 

The paper is organized as follows. In Sec.|TT]the Maxwell equations for two coupled waveguides system are solved, 
and, in particular, expressions are obtained for the fields and the dispersion relations. In Sec. IHII the Maxwell stress 
tensor technique is used to calculate the optical force, whose behavior is studied in detail as a function of the distance 
between the waveguides and of the incident light frequency. In Sec. IIVI we quantitatively estimate the mechanical 
deflection produced by the radiation force in a realistic air-bridge device and we show that it is strong enough to 
be experimentally measured by means of Atomic Force Microscopy techniques. Finally, in Sec. |V] we derive closed 
expressions for the force in the limiting cases of small and large separation. 



II. GUIDED ELECTROMAGNETIC MODES IN COUPLED SLAB WAVEGUIDES 



The physical configuration we consider in the present paper is shown in Fig[TJ a pair of parallel dielectric slabs of 
thickness s and refractive index rig separated by a distance 2a and embedded in a host medium of lower refractive 
index nu < ng. The axis x is orthogonal to the slabs and light is assumed (without loss of generality) to be in a plane 
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FIG. 1: Left panel: Schematic view of the two coupled waveguide system under study. The red arrows represent the input laser light. 
Right panel: Transverse cut of the system along a fixed-2 plane. The TE/TM arrows indicate the polarization of the electric field in the 
two polarization states. 



wave state propagating along z. Under this assumption, the electric and magnetic fields can be written in the form: 



E(r,t) = Re E(x)e 



i(ut—(3z) 



H(r, t) = Re 



H(x) e 



i(ujt—(3z) 



(1) 
(2) 



[3 being the wave number of the propagation along z and to the angular frequency. The y dependence of the fields has 
disappeared thanks to the translational invariance of the system in the waveguide plane, and to the choice made for the 
direction of propagation. The general solution for the electromagnetic fields can be obtained as a linear combination 
of such plane waves. 

Depending on the orientation of the fields with respect to the propagation direction, two independent electromag- 
netic polarizations states can be identified, known as Transverse Electric (TE) and Transverse Magnetic (TM) 0, [y] . 
The TE polarization state is characterized by E x = E z — H y = 0, while the TM polarization is characterized by 
E y — H x = H z = 

For each TE/TM polarization state, the field wavefunction is determined by solving the corresponding Maxwell 
equations. Since we are considering modes which are guided by the waveguide system, the electromagnetic field is 
confined in the slabs, where the wave- vector along the x-direction is purely real k x = ±k 



±Vk 



2 ?i 2 



/3 2 , with 



k = bj/c. In the surrounding host medium the electromagnetic field is evanescent with a purely imaginary wave- vector 
K = ±ia = ±iy/f3 2 - k 2 



H' 



For the TE/TM polarizations, this is summarized in a field wavefunction (that is the E y or H y fields for the 
TE/TM modes, respectively) which reads: 



Ey , Hy — A < 



a± e 



COs[k(x — a) + <f)2\ 

a?2 cosh(ax) + a' 2 smh(ax) 
cos[k(:e + a) — fa] 
a 3 e +a( - x+a+ ^ 



(x > a + s) 
a < x < a + s 
—a < x < a 

—a < x < —a 
x < —a — s, 



(3) 



The system having reflection symmetry with respect to the x = plane, the electromagnetic modes can be classified 
as symmetric and antisymmetric depending on the symmetry operation of the field wavefunction, namely E y for the 
TE modes or H y for the TM ones. Note that this classification is related, but not identical to the usual one in terms 
of the reflection symmetry of the full electromagnetic field, where the magnetic H field transforms as a pseudovector, 
while E is instead a vector. 
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The dispersion relation connecting (3 to to, as well as a relation between the amplitude coefficients a and the phases 
cj) in the different regions are then obtained by matching the field §3§ in the different regions according to the symmetry 
of the electric and magnetic fields. For the symmetric modes, the dispersion law has the following analytical, yet 
implicit form: 



ks = arctan 
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(4) 
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for the TE/TM polarizations, respectively. The dispersion laws for the anti-symmetric modes are obtained from 
the symmetric ones by replacing tanh(x) i— > 1/ tanh(a;). In the following, we shall see that this fact holds for other 
physical quantities as well. In all these dispersion laws, the (integer) quantum number m > specifies the number of 
nodes of the wavefunction inside each slab. 

In summary, the optical modes propagating along z in a given two waveguide system are classified by their TE/TM 
polarization state, the symmetric/anti-symmetric character of the field wavefunction with respect to reflections on 
the x = plane, and the number m > of nodes inside each waveguide. Thanks to the scaling properties of the 
Maxwell equations, the qualitative shape of the dispersion relations depends on the geometrical parameter a/s only, 
which quantifies the ratio between the spacing a and the thickness s of each slab. The absolute value of s fixes the 
natural frequency scale lo s — 2?t c/s. 

In Fig. [2j we have plotted the dispersion of the different modes (TE/TM, symmetric and anti-symmetric) of the 
coupled slab waveguide system for different values of the geometrical parameter a/s. All the branches have a lower 
cut-off to the frequencies that can actually propagate in a given mode of the coupled waveguide system. At the cut-off 
point, the waveguide dispersion coincides with the free photon dispersion oj = c(3/uh of the host medium. 

At infinite separation a/s — oo, the dispersion reduces to the one of an isolated waveguide, and for every polarization 
state the symmetric and anti-symmetric branches are degenerate. As the two waveguides are pushed closer, this 
degeneracy is lifted, and every branch experiences a frequency shift whose sign depends on its symmetric or anti- 
symmetric nature. As usual for the bonding/anti-bonding electronic states in diatomic molecules [7|, the symmetric 
states are pushed toward lower frequencies by the coupling, while the anti-symmetric ones are pushed toward higher 
frequencies. As a consequence, the cut-off frequency experiences itself a shift of the same sign and comparable 
magnitude. 

In the next sections, we shall study the electromagnetic pressure acting on each of the two slab waveguides because 
of the presence of the other slab. The pressure being proportional to the intensity of light propagating along the 
waveguide system, it is important to relate the amplitude coefficient A in (J3]) to the power density P for unit length 
in the y-direction. This is easily calculated from the flux of the Poynting vector through a planar section orthogonal 
to the propagation direction. 

For the symmetric modes, we obtain 
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for the TE/TM polarizations, respectively. The expressions for the corresponding antisymmetric modes are obtained 
by replacing tanh(x) i— > tanh(cc) -1 in ^ and ([7]). 



III. FORCE BETWEEN TWO PLANAR WAVEGUIDES 



Starting from the electromagnetic field profiles discussed in the previous section, we shall now proceed with a 
calculation of the average electromagnetic force F acting on the slabs when a monochromatic wave of frequency lu is 
propagating along the coupled waveguide system in a well-defined mode. To keep the treatment as simple as possible, 
we shall assume this force to be equilibrated by some other, unspecified, force which keeps the system at mechanical 
equilibrium at all times. The calculation of the force will then be performed in the framework of the macroscopic 
electrodynamics of continuous media using the Maxwell stress tensor T [H, Q . This technique allows one to directly 
calculate the force, and has been extensively used in the literature to estimate forces of electromagnetic origin in 



4 



High energy limit (well guided modes) 



0.2840- 




0.2 0.4 0.6 0.8 1.0 



o 
c 

CD 
Z3 
O" 



(3 [2tt/s] 



0.05 



0.10 



0.15 



0.20 



0.25 



P [27I/S] 



FIG. 2: Main panel: dispersion relation of the different symmetric and antisymmetric branches of the fundamental (m = 0) and first 
order (m = 1) TM mode for various waveguides separations. The thick continuous lines refer to the infinite separation a/s — > co case 
where symmetric and anti-symmetric modes are degenerate. Upper and lower lines with respect to the infinite separation one refer to 
symmetric (Sym) and antisymmetric (Asym) modes respectively; the thin continuous lines are for a/s = 0.5, while the dashed ones are 
for a shorter separation a/s = 0.25. In the smaller panels on the right we have highlighted the behavior of the m = modes, respectively 
far from the cut-off (high energy limit) and near the cut-off (low energy limit). 



many contexts, from clusters of dielectric spheres [T(| to waveguides 0], to quantum fluctuations as in the Casimir 
effect [ll] ) . An important point of this formalism is that it does not make any assumption on the microscopic nature of 
the material media under examination and therefore can easily take into account absorption effects. In the following 
we shall focus our attention on the case of dielectric slabs with a real refraction index ns = n embedded in air, for 
which Tin — 1- 

As all the fields have a monochromatic time dependence at frequency to, momentum conservation arguments show 
that the average electromagnetic force acting on a body is equal to the surface integral of the time averaged Maxwell 
tensor 



T, 



Re 



E. L E* + H t H* - - Sij i E k E t + B k B* k 



(8) 



over an arbitrarily chosen closed surface S enclosing the body with outward orientation: 



fij(r)da 



3' 



(9) 



In our specific configuration, a good choice for £ is the one shown in FigQJ), that is a cylinder of axis parallel to z 
and with rectangular bases parallel to the xy plane. 
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FIG. 3: Electromagnetic pressure for the fundamental (m = 0) TE (continuous lines) and TM (dashed lines) modes for s = 310 nm thick 
silicon (n = 3.5) waveguides, a wavelength A = 1.55/im and a power density of P = 20mW//im. Upper (lower) curves, corresponding to 
a repulsive (attractive) force, refer to antisymmetric (symmetric) modes. 



Thanks to the reflection symmetry of the whole set-up with respect to the xz plane, the y component to the 
force vanishes F y — 0. Also along the light propagation direction z the force is zero, indeed as the geometry of the 
waveguide system is symmetric with respect to the xy plane, and the dielectric medium is non-absorbing Im[n] = 0, 
the combination of this reflection symmetry and the time-reversal is a symmetry of the problem. The electromagnetic 
force is therefore directed along the x direction. The contribution of the two planar sides parallel to the xz plane 
cancel each other by translational symmetry, as well as the one of the two bases parallel to the xy plane. We are 
therefore left with the sides parallel to the yz plane; the x component of the force due to these sides only involves 
the xx component of the Maxwell stress tensor; because of the translational symmetry of the configuration under 
examination, this quantity can only depend on the x coordinate: 

f xx = -| [|£y 2 + 1£ 2 | 2 - 1^| 2 + c^i (|^| 2 + 1^| 2 - i^i 2 )] . (io) 
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Inserting the explicit form of the fields, it is immediate to see that T xx = in the region \x\ > a + s external to 
the waveguide system. This is due to the evanescent wave character of the field in this region. The electromagnetic 
pressure p (i.e., the force per unit length along z and unit width along y) is therefore equal to —T xx evaluated in the 
region between the two waveguides, |z| < a. Positive (negative) signs for p respectively indicate repulsive (attractive) 
forces between the waveguides. Plugging in (|10j) the explicit expression of the fields found in the previous section 
(Sec. HII, we obtain the following results for the symmetric TE/TM modes: 



The expression for the TE/TM antisymmetric modes are again found by replacing tanh(x) i— > tanh(a;) _1 in (|1 1|) 
and (|12[) . In the following, we will work at constant laser frequency lj and power density P, so that the amplitude 
coefficients Ate and Atm have to be obtained by inverting ([6]) and ([7]). 

Note that since the Maxwell stress tensor is bilinear in the local fields, the effect described here does not rely on the 
coherence of the light beam, thus the results of the present paper directly extend to the case of an incoherent, thermal 
source. Indeed the total pressure induced by a source with spectral distribution f(u>) is simply J dui f(ui)p(u>), where 
p is given by Eqs. CU} and lfl2]). 

Regarding the monochromatic source, considered hereafter, one has to distinguish two cases, depending on whether 
the laser frequency is far from or close to the cut-off frequency of the considered mode. Results for the first case are 
shown in Figure [3] for the m = mode in 310 nm thick silicon [n = 3.5) waveguides: the pressure is plotted as a 
function of the separation 2a between the waveguides. Fixed values are taken for the power density P — 20mW//Ltm 
and the wavelength A = 1.55 /im of the wave. The main feature is that the pressure is always attractive for both 
TE (continuous lines) and TM (dashed lines) symmetric modes, while it is repulsive for the spatially antisymmetric 
modes. Formally, this is an immediate consequence of (fTTj) and (fT2"| . once one takes into account the fact that for 
guided modes > k. In magnitude, the force is a monotonically decreasing function of the separation a. Physically, 
this behavior has an immediate explanation in terms of the analogy with two coupled well models: as shown in FigJ^l 
the frequency of the antisymmetric (symmetric) mode for a given (3 monotonically grows (decreases) as the waveguides 
are brought closer. 

At large distances, the decay of the pressure as a function of distance is exponential and for a given TE/TM 
polarization, the symmetric/anti-symmetric modes only differ by the sign of the pressure. Since the TE mode is more 
confined in the slabs than the TM one, it has a shorter characteristic length of the exponential. At short distances, it 
is interesting to note that the TM symmetric mode produces an significantly enhanced pressure as compared to the 
corresponding TE one. A physical explanation of this behavior is readily obtained by comparing the E x {x) electric 
field profiles of the m = symmetric TE and TM modes, as shown in FigfJ] while the TE mode profile has a smooth 
spatial dependence, the TM one is strongly concentrated in the region between the two slabs. This feature is typical 
of TM modes [l2[ , and originates from the continuity of the normal component D z of the electric displacement vector 
at the slab interface, which introduces a n 2 factor between the values of E z at the two sides of the interface. 

If the laser frequency is not far from the cut-off of the mode, the dependence of the pressure p on the separation a 
is somehow richer. In Fig[S^ we consider the case of a thicker waveguide s sa I fim. For the wavelength A = 1.55 /im 
under consideration, all the modes up to m = 3 are well confined, while we are just above the isolated (a/s = oo) 
waveguide cut-off for the m = 4 mode. Since the cut-off frequency of anti-symmetric modes increases for decreasing 
a, it exists a cut-off separation a co below which light of the given wavelength ceases being guided in the m — 4 mode. 
When a — > a co from above, the spatial size cr _1 of the mode diverges, and the field amplitude between the guides 
tends (for a given power density P) to zero. As a consequence, the pressure p initially grows for decreasing a, attains a 
maximum value at some separation value, and then goes back to zero as the cut-off separation a co is approached [l3|. 
Clearly, the cut-off separation is larger for thinner waveguides (Fig[5jD). This behavior can also be explained in terms 
of the two coupled well model: when the eigenstates of the independent wells are close to the continuum threshold, 
there exists a value of the distance (i.e., of the coupling strength), at which the antisymmetric state ceases to be bound 
and enters the continuum. On the other hand, for the tightly confined m < 3 modes, the physics is qualitatively the 
same as in FigJ3J 




TE 



(11) 




TM 



(12) 



7 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 

Position [|im] 

FIG. 4: Profile of the E x component of the electric field for the symmetric TE and TM m = modes. 




FIG. 5: Plot of the pressure as a function of the separation a for (a) modes with a different order m; (b) for the same m = 4 and slightly 
different waveguide thickness. Upper (lower) curves, corresponding to a repulsive (attractive) force, refer to antisymmetric (symmetric) 
modes. 



IV. RADIATION INDUCED DEFLECTION: THE CASE OF THE AIR-BRIDGE WAVEGUIDE 

Although the value of the radiation pressure obtained in the previous section might seem at a first glance rather 
small, it can have a observable effect in nanometric devices. In order to provide a quantitative estimate of such effects 
we consider the mechanical deflection induced by the optical force on an air-bridge double slab waveguide system 
made of two thin silicon slabs of length L whose opposite edges are clamped to the substrate. Such a device can be 
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FIG. 6: Maximal value £max of the deflection £ produced by the TM symmetric mode (with power P = 20mW//jm and wavelength 
A = 1.55/im) of the radiation induced force on an silicon double air-bridge waveguide of thickness s = 310nm, as a function of the initial 
separation 2a. Three different values of the slabs length are considered: L = 30/^m (solid), L = 35^tm (dashed) and L = 40/mi (dotted). 



realized, e.g., via chemical etching technique in [l5j . 

The deflection £ induced by the radiation pressure on the device can be evaluated using the Euler-Bernoulli beam 
equation [l4j 

o4 p 

EI d^ = MC) ' (13) 

where, £ is the separation between the slabs, E is the Young modulus, and / = (l/12)/is 3 is the area moment of 
inertia of the slab's cross section, whose width and thickness are h and s. The maximum deflection ^ max induced on 
each silicon {E — I69GPa) slab by a power P = 20mW//jm at a wavelength A = 1.55^m propagating in the TM 
symmetric mode is shown in Fig.® as a function of the initial separation 2a for three different values of the slab 
length L = 30, 35, 40 fim, and for slab's thickness s — 3I0nm. For this parameter choice, a linearized treatment of 
Eq. tfTB")) around the initial slab distance already provides accurate results. 

From figure © one can see that the maximal value £ max of the deflection £ induced by the optical force can reach 
the values of tens nanometers for reasonable structural parameters: slabs lengths 30/im< L < 40/xm (the deflection 
is strongly sensitive to L), and waveguides separations 50nm< 2a < 200nm. Such deflections are of the same order of 
previous investigations made for different configurations [3| and are accessible to standard Atomic Force Microscope 
techniques [lq |. As a final remark, as the pressure is roughly speaking inversely proportional p cx l/v g to the group 
velocity v g of the mode, one expects that the radiation pressure, and hence the mechanical deflection, can be enhanced 
if the v g is slowed down by a longitudinal modulation of the structure as in coupled resonator optical waveguides [13] • 



V. LARGE AND SMALL DISTANCE BEHAVIOR 



In this section we derive the asymptotic behavior of the pressure for small and large slab separation, which provides 
a deeper insight on the findings of the previous section. 
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A. Large distance behavior 

For large separation a/s —* oo, the force has the typical exponential decay of two-well systems in the tight-binding 
limit. As long as the modes are confined, the general qualitative trend is that the larger the characteristic length 
(o'oo E/ ' TM ) _1 , the larger the field in between the slabs and consequently the stronger the force. More specifically: for 
a given order m the pressure is (slightly) stronger for the TM mode than for the corresponding TE one. (ii) The 
pressure is stronger for higher m modes. 

Quantitative expressions for the case of symmetric (anti-symmetric) modes can be obtained by expanding Eqs. 
and (fT2"j) for large distances: 



TM,± , tj 
Ptb = ±P 
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fc 2 



c(3, 



TM, 



-2aoi 



(14) 



(15) 



, TE/TM TE/TM a TE ITM , , , . c .. .. , 

where Koo , Uoa , Poo are evaluated at infinite separation a/s 
the symmetric (antisymmetric) mode. 



oo. As before, the +(— ) sign refers to 



B. Small distance behavior 



In order to get some analytical insight in the small distance regime, it is useful to expand all the waveguide 
parameters in powers of the slab separation a, while keeping P and lo constant. Let us consider the specific case of 
/3(a) = /3q + Pi a + 0{a 2 ). The zeroth order /3q = (3(a = 0) is the wave number of the propagation along a single 
waveguide of double thickness 2s, and has to be calculated from the dispersion laws once the hyperbolic tangent is 
replaced by its limiting value 1 and the doubled thickness is taken into account as s — > 2s. The first order term /3i is 
given, e.g., for m — symmetric modes by: 



pi(a) = - 
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(16) 
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Note that (3i is negative for all modes, indicating that the force is attractive. Along the same lines, analytical 
expansions can be obtained for the pressure p. The a = value p for the symmetric TE/TM modes has the form 
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Analogous algebra leads to the corresponding expressions for the anti-symmetric TE/TM modes, which have the 
form: 
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Starting from these formulas, a physical explanation can be provided for the remarkable facts observed in Fig[2] for 
the m = modes, namely the suppressed value of the force for the anti-symmetric TM mode and the enhanced value 
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for the symmetric TM mode with respect to the TE modes. As long as we are considering well confined modes, 
the (To's of all modes have almost comparable values, somewhat larger than the Kq's. This explains the general fact 
that the force is about a factor 2 weaker for the TE antisymmetric mode than for the corresponding symmetric one. 
The behavior for the TM modes can be explained starting from the value n — 3.5 of the refractive index, which is 
significantly larger than 1: thanks to the n i in the denominator, the force p M '~ for the antisymmetric TM mode is 
dramatically suppressed of a factor rj 150 (for the chosen value n — 3.5 of the refractive index) with respect to the 
one p ' + for the symmetric TM mode. For similar reasons, the force p ' + is enhanced of a factor m 13 over p ,+ 
because of the n's in the denominator (see Fig. [3]). 



VI. CONCLUSIONS 

In conclusion, we have performed an analytical study of the optical forces appearing between a pair of parallel slab 
waveguides when light is propagating through them. Depending on the spatial symmetry of the mode wavefunction, 
the sign of the force can be either attractive or repulsive. The dependence of the force as a function of the separation 
between the slabs has been characterized for the different polarization states, and analytical expressions have been 
obtained for both the large and the small distance limits. A strong enhancement of the force has been identified for the 
symmetric TM mode, as well as a suppression for the antisymmetric TM one. Simple physical explanations have been 
provided for these features. A quantitative study for typical air-bridge configurations confirms that the mechanical 
deflection of the structure induced by the optical force can be measured by standard Atomic Force Microscopy 
techniques. 
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